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INTRODUCTION 

Quantum spin systems show a very rich variety of fas- 
cinating phenomena. Despite the difficulty of describing 
such systems due to the large number of degrees of free- 
dom, the development of both analytical and numerical 
techniques have allowed us to investigate and understand 
many of those phenomena. This success story is due, in 
part, to the discovery of renormalization group methods 

m which are very powerful scalable numerical tech- 
niques which seem to describe amazingly well the ground 
states of ID spin systems with local interactions in terms 
of matrix product states (MPS) HQ. 

On the other hand, one may also investigate many- 
body quantum systems with a quantum computer or a 
quantum simulator la]. For instance, one may use an 
adiabatic algorithm j6| to prepare the quantum simula- 
tor in the ground state of a prescribed Hamiltonian Q, 
and then measure all one- and two-particle correlation 
functions. It seems that this approach may allow us to 
investigate the properties of many-body quantum sys- 
tems that otherwise would not be describable by ordi- 
nary computers using existing numerical algorithms. In 
fact, an important effort is being made in order to build 
quantum simulators for some specific tasks ||. 

In this paper we present a series of analytical results 
which quantify the accuracy of the renormalization group 
methods. First, we give an expression that bounds the 
errors made by approximating general states of N spins 
by MPS in terms of the Renyi entropies corresponding to 
different subsystems. This result, when combined with 
the scaling of those Renyi entropies for critical ID sys- 
tems, implies that both short- and long-range properties 
of the low energy states of critical spin chains can be very 
precisely described in terms of MPS with a number of pa- 
rameters which only scales polynomically with N. This 
indicates that for those systems a quantum computer or 
simulator may not give us a big advantage with respect 
to a classical one. Second, we determine a bound on the 
error made in the approximation of 2-particle reduced 
density operators coming from translationally invariant 



states in ID by MPS in the limit N — > oo. Finally, 
we illustrate the accuracy of renormalization group algo- 
rithms by some examples. 



GROUND STATES AS CONVEX PROBLEMS 

The determination of the ground state of a quan- 
tum system consisting of interacting spins on a lattice 
is highly nontrivial due to frustration effects. If one con- 
siders e.g. a Heisenberg antiferromagnetic interaction 
between neighboring spins on a 1-D spin chain, the en- 
ergy would be minimized if all the reduced density opera- 
tors of all neighboring spins would correspond to singlets. 
The monogamy properties of entanglement [!J, [lfj, [Tl], Il3 | 
however imply frustration effects which forbid the exis- 
tence of such a state: if a particle A is maximally en- 
tangled with B, then it cannot be entangled with C. 
From the technical point of view, this can readily be 
proven by showing that the convex set of states obey- 
ing the semidefinite constraints Pabc > 0;Trc (pabc) = 
|S)(S| =Tr B (pabc) with \S) = (|01)-|10))/V5 is empty. 
The frustration effects become stronger when the coor- 
dination number of the lattice increases 0, 0] . By in- 
voking the quantum de-Finetti theorem one can prove 
that in the limit of large coordination number all possi- 
ble ground states will be separable 0, 0] , which on its 
turn implies that mean-field theory becomes exact in the 
limit of infinite dimensional lattices. 

Due to the variational nature of ground states, there 
always exists a ground state with the same symmetries 
as the associated Hamiltonian. If the Hamiltonian has 
translational symmetry and consists of 2-body nearest 
neighbor interactions, then it is clear that the energy 
of a state with the right symmetry is completely deter- 
mined by its reduced density operator of 2 neighboring 
spins. The reduced density operators arising from these 
(eventually mixed) states with a given symmetry form a 
convex set, and the energy for a given Hamiltonian will be 
minimized for a state whose reduced density operator is 
an extreme point in this set. More specifically, the equa- 
tion Tr(Hp) = E determines a hyperplane in the space of 
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reduced density operators of states with a given symme- 
try, and the energy will be extremal when the hyperplane 
is tangent to the convex set (E = E extr ). The problem 
of finding the ground state energy of nearest neighbor 
translational invariant Hamiltonians is therefore equiv- 
alent to the determination of the convex set of 2-body 
reduced density operators arising from states with the 
right symmetry. Strictly speaking, these two problems 
are dual to each other. In the case of quadratic Hamilto- 
nians involving continuous variables, the determination 
of this convex set was solved for fairly general settings in 
[l4| by means of Gaussian states. The determination of 
this convex set in the case of spin systems however turns 
out to be much more challenging. 

Let us illustrate this with a simple example. Consider 
the XXZ-Hamiltonian with nearest neighbor interactions 

n = - sfsj + sfsy + as?s? 

<i,j> 

on a lattice of arbitrary geometry and dimension [36j|. 
Due to the symmetries, the reduced density operator of 
two nearest neighbors can be parameterized by only two 
parameters [23 : 

p = ~ (1 ® 1 + x(a x <g> a x + a y <E> <J y ) + za z <g> a z ) . 

Positivity of p enforces —\<z<l — 2\x\, and the state 
is separable iff 1 + z < 2\x\. In the case of an infinite 
1-D spin chain, the ground state energy E{A) has been 
calculated exactly and this determines the tangent 
hyperplanes 

2x + zA + E(A) = 

whose envelope makes up the extreme points of the con- 
vex set of reduced density operators of translationally 
invariant 1-D states: the boundary of this convex set is 
parameterized by 

z = -dE(A)/dA 

x = -(E{A) + dE(A)/dA)/2, 

which we plotted in Figure 1. We also plot the boundary 
for the 2-dimensional square lattice. These 2-D data were 
obtained by numerical methods 0, 0, 0|); of course 
this convex set is contained in the previous one, as all 
the semidefinite constraints defining the set correspond- 
ing to 1-D are strictly included in the set of constraints 
for the 2-D case. Finally, we plot the set of separable 
states, which contains the reduced density operators of 
the allowed states for a lattice with infinite coordination 
number. The boundary of this separable set is given by 
the inner diamond; this immediately implies that the dif- 
ference between the exact energy and the one obtained 
by mean field theory will be maximized whenever the hy- 
perplane that forms the boundary of the first set will be 




Figure 1: Convex sets of the possible reduced density opera- 
tors of translational invariant states in the XX-ZZ plane: the 
big triangle represents all positive density operators; the in- 
ner parallellogram represents the separable states; the union 
of the separable cone and the convex hull of the full curved 
line is the complete convex set in the case of a 1-D geometry, 
and the dashed lines represent extreme points in the 2-D case 
of a square lattice. The singlet corresponds to the point with 
coordinates (-1,-1). 



parallel to this line. This happens when A = — 1 (in- 
dependent of the dimension!), which corresponds to the 
antiferromagnetic case, and this proves that the "entan- 
glement gap" 0] m the XXZ-plane is maximized by the 
antiferromagnetic ground state for any dimension and 
geometry. Similarly, it proves that the ground state is 
separable whenever A > 1 and A = — 00. Note also that 
in the 2-D case, part of the boundary of the convex set 
consists of a plane parameterized by 2x + z + E(l) = 0. 
This indicates a degeneracy of the ground state around 
the antiferromagnetic point, and indicates that a phase 
transition is occurring at that point (more specifically 
between an Ising and a Berezinskii- Kosterlitz-Thouless 
phase). It would be very interesting to investigate this 
further, and it is fascinating that this phase transition 
could be detected by just looking at the structure of these 
low-dimensional convex sets. 

The previous discussion implies that ground states of 
quantum spin systems are very special: they are com- 
pletely determined by their 2-body reduced density op- 
erators. Typically, it even holds that the ground state 
of an interacting spin system is unique, which implies 
that most extreme points in the convex set of reduced 
density operators uniquely characterize a state with the 
right symmetry properties. This is very good news if 
we want to create families of variational ground states: 
it suffices to approximate well the local properties of all 
translational invariant states. The family of matrix prod- 
uct states (MPS) 0,0, H(J and generalizations to higher 
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Figure 2: A one dimensional spin chain with finite correlation 
length ^ C orr\ Iab denotes the distance between the block A 
(left) and B (right). Because Iab is much larger than the 
correlation length £ C orr, the state pab is essentially a product 
state. 

dimensions (PEPS) 0,0, H3 were exactly created with 
this property in mind; the amazing accuracy of renor- 
malization group algorithms is precisely related to the 
fact that the convex set under consideration can be very 
well approximated with the reduced density operators of 
MPS. Both Wilson's numerical renormalization group Q 
and DMRG [2J, |23| can indeed be reformulated as varia- 
tional methods within the MPS [iEof. 

The reason why ground states of gapped quantum spin 
systems can be parameterized as MPS/PEPS can be un- 
derstood by the following handwaving argument. A cele- 
brated theorem of Hastings 01 states that in the case of 
a gapped system all correlations are decaying exponen- 
tially. Let us therefore consider a 1-D gapped quantum 
spin system with correlation length £ corr . Due to the fi- 
nite correlation length, it is expected that the reduced 
density operator pab obtained when tracing out a block 
C of length Iab ^S> £ C orr (see figure 01 is equal to 

Pab ^ Pa® Pb (1) 

up to exponentially small corrections. As noted by T. 
Osborne, Hastings theorem does not imply the validity 
of equation , as it was shown in [2j| that orthogonal 
states exist whose correlation functions are exponentially 
close to each other; although it would be very surprising 
that ground states would exhibit that property, this pro- 
hibits to turn the presented handwaving argument into a 
rigorous one. But let's for the sake of the argument as- 
sume that equation Q is true. Then, the original ground 
state \ipABc) is a purification of this mixed state, but it 
is of course also possible to find a purification of the form 
iV'ACi) ® \ipBC r ) ( U P to exponentially small corrections) 
with no correlations whatsoever between A and B; here 
Ci and C r together span the original block C. It is how- 
ever well known that all possible purifications of a mixed 
state are equivalent to each other up to local unitaries on 
the ancillary Hilbert space. This automatically implies 
that there exists a unitary operation Uq on the block C 
(see figure |2J that completely disentangles the left from 
the right part: 

t A ®U C ® &c\tpABC - \i>AC,) ® \lpBC r )- 
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Figure 3: Convex sets in the XXZ-plane: the inner diamond 
borders the set of separable states (see Fig. 1). Dash-dotted: 
extreme points of the convex set produced by MPS of D = 2. 



This implies that there exists a tensor A l a a with indices 
1 < a, (3,i < D (where D is the dimension of the Hilbert 
space of C) and states IV'-F)' 1^9 ) defined on the 

Hilbert spaces belonging to A,B,C such that 

Wabc)^J2 A *MM?M)- 

Applying this argument recursively leads to a MPS (to 
be defined in equation an d indeed gives a strong 
hint that ground states of gapped Hamiltonians are well 
represented by MPS. In this paper we will prove that this 
is indeed true, even in the case of critical systems. 

As a good illustration of the actual accuracy obtained 
with MPS, we calculated the convex set obtained with 
MPS in the thermodynamic limit for the XXZ-chain with 
D = 2, where D is the dimension of the matrices in the 
MPS. It is almost unbelievable how good the exact con- 
vex set can be approximated. Note that typical DMRG 
calculations have D ~ 200, and that the accuracy grows 
superpolynomial in D. Note also that the D = 1 case 
corresponds to mean-field theory, whose corresponding 
convex set coincides with the set of separable states. 

The same argument involving the notion of a cor- 
relation length applies in higher dimensions and indi- 
cates that PEPS represent ground states of gapped local 
Hamiltonians well. Note however that the convex set in 
the 2-D case is much closer to the separable one than 
in the 1-D case; this gives a hint that PEPS of smaller 
dimension will suffice to give the same accuracy as in 
the 1-D case. In the next section we will quantitatively 
bound how well a translationally invariant state can be 
represented in terms of a MPS, and will analyze the corre- 
sponding implications for the description of ground states 
of ID spin chains. 
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ANALYTICAL BOUNDS FOR MATRIX 
PRODUCT STATES 



Let us start by recalling the expression of D- 
dimcnsional MPS describing N spins of dimension d, 



Proof: We can always write ip as a MPS of dimension 
D = 2 N / 2 and fulfilling J3J. Let us now consider the D- 
dimensional MPS ^d) which is defined by the D x D 
matrices 

.d,i...d) (i- e - the upper-left block of 
A^ 1 ). The goal is now to bound (tp\tpr>). The gauge 
conditions were chosen such as to make the task simple: 



\4>d) 



E 



tr 



\h,...,i N ). (2) (V> D |V> =Tr[$ 2 (---$ 



»1 tN = l 



Here, A^" 1 , . . . , A^' d are, in general, Dk x Dk complex 
matrices, with Dk+i = Dk < D. 

We will show first how well one can describe general 
states in terms of MPS. We will choose D\ = Dn = 1, 
i.e., A^ n and A^ m will be vectors. In this case, one 
can impose the following gauge condition |2ri |2?1 | : 

^^[m]j^[m]if _ ^ ^[ro]if A [ro+l]^[m]i = \l m ] , 

i i 

(3) 

Here Al m ] represents the diagonal matrix with the corre- 
sponding eigenvalues {Al™! 2 } sorted in decreasing order, 
i.e. Al™! 1 > A'™' 2 > . . .. We consider an arbitrary state 
l^} and denote by 

{ M Hi }, l = i..Af a = r in(a ' JV " a) 

the eigenvalues of the reduced density operators 

p a = Tr a+ i, Q , +2 ,...,N|V , )(V>l> 

also sorted in decreasing order. 

In the following we will investigate: (i) How well a 
general quantum state can be approximated by a D— 
dimensional MPS; (ii) How well the reduced density op- 
erators of a translationally invariant state can be describe 
by the one corresponding to a Z?-dimensional MPS. 



General states: analytic results and implications 

We start out with the main result of this paper, which 
gives an upper bound to the error made by approximating 
a general state by a MPS. As we will show below, this 
has very important implications in the performance of 
the renormalization algorithms to describe ground states 
of ID spin chains. 

Lemma 1 There exists a MPS \iPd) of dimension D 
such that 



N-l 



\H)-Hd)\\ 2 <2Y,^(D) 



a=l 



jv-2 (%n-i (A N ~ 1 P) P)P---)P]; 

(4) 

here P = Ek=i\ k )( k \ and $ m {X) = ^^H't^H- 
represents a trace-preserving completely positive map 
(TPCP-map) parameterized by the Kraus operators 
^[ro]i_ us now recursively define 



YW = $, 



_ A [iv-i] p. 



observe that A^J = $ fe (Al fe+1 l). We want a bound on 
Tr|AW - Yi|, as equation © is equal to Tr(Y^). The 
crucial property we need is that TPCP-maps are con- 
tractive with relation to the trace-norm [3^|: Tr|$(X)| < 
TrlXl. It follows that 



Tr | A H _ Y [k]| _ Tr | $k / A [k+i] _ Y [ k +i]p 

< Tr | A [k+i] _ Y [ k+1 lp| 

< Tr|A[ k+1 l - Yl k+1 l | + Tr|A[ k+1 l (1 - P)| 



< 



Note that the last term in the sum is exactly given by 

Y^a=D+i A[ fe+1 ' Q . The theorem now follows immediately 
by recursion and by observing that (iPdI^d) < 1 by sim- 
ilar arguments. □ 

The implications of this result are very strong: it 
shows that for systems for which the e a (D) decay fast 
in D, there exist MPS with small D which will not only 
reproduce well the local correlations (such as energy) 
but also all the nonlocal properties (such as correlation 
length). The following lemma now relates the derived 
bound to the Renyi entropies of the reduced density op- 
erators, through which one can make the connection to 
the ground states of ID Hamiltonians. The Renyi en- 
tropies of p are defined as 



s a ( P ) 



1 



1 



log (Trp") 



where e a (D) = £iL°b+i p [a 



and we will consider < a < 1. We denote as before 
e(D) = w ith Xi the nonincreasingly ordered 

eigenvalues of p. Then we have 

Lemma 2 Given a density operator p. If < a < 1, 
then log(e(D)) < (s a (p) - log ^) . 

Proof: Let us first characterize the probability distri- 
bution that has maximal possible weight in its tail (i.e. 
P = ESd+iP«) f° r a gi ven Renyi-entropy. Introducing 
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a free parameter < h < (1 —p)/D, such a probability 
distribution must be of the form 

Pi = l-p-(D-l)h 

h = P2=P3 = ■■ -PD+p/h 
PD+p/h+li ' "Poo = 

because this distribution majorizes all other ones with 
given p,D,pd (Renyi-entropies are Schur-convex func- 
tions). For a given p, D, h, it holds that 

J2p? = (l- P -(D-l)h) a + (D-l+p/h)h a 

i 

> DhT+ph - 1 . 
Minimizing this expression with relation to h, we get 

52tf>(D 1 - a p a )/((l-a) 1 - a a a ). 

i 

Denoting S a (p, D) the minimal possible entropy for given 
p, D, we get 

1 / D 1 ~ a n a \ 

S«( P ,D)>— ■° 8 ( (1 _ tt) ,! aa „ ) 

and hence 

P <exp(i^(^)-lo g ^)). 

The proof now follows by replacing S a (p,D) by S a (p). 

□ 

This lemma is very interesting in the light of the fact 
that in the case of critical systems, arguable the hardest 
ones to simulate [j^J, the Renyi-entropy of a contiguous 
block of L spins scales as [H H El H2 

S"(p £ )~£^fl + lW(L) (5) 

for all a > 0; here c is the central charge. The fact 
that the eigenvalues of pL decay fast has previously been 
identified as a indication for the validity of the DMRG- 
approach [2^ . The truncation error [3, , which has 
been used in the DMRG community as a check for con- 
vergence, is essentially given by e(D) — e(2D) and there- 
fore indeed gives a good idea of the error in a simulation. 

Let us investigate how the computational effort to sim- 
ulate such critical systems scales as a function of the 
length N — 2L of the chain. Let us therefore consider the 
Hamiltonian associated to a critical system, but restrict 
it to 2L sites. The entropy of a half chain (we consider the 
ground state \ip ex ) of the finite system) will typically scale 
as in eq. 10 but with an extra term that scales like 1/N. 
Suppose we want to enforce that || \ip ex ) — \i/jd) || 2 < eo/L 
with eo independent of L |4d| . Denote the minimal D 
needed to get this precision for a chain of length 2L by 



Dl. Following lemma (1) and the fact that the entropy 
of all possible contiguous blocks reaches its maximum in 
the middle of the chain (hence p < £q/L 2 is certainly 
sufficient), lemma (1) and (2) combined yield 

D L <cst[- — L^-. 

V(l - a)e J 

This shows that D only has to scale polynomially in L to 
keep the accuracy eo/L fixed; in other words, there exists 
an efficient scalable representation for ground states of 
critical systems (and hence also of noncritical systems) 
in terms of MPS! Such a strong result could not have 
been anticipated from just doing simulations. 

Now what about the complexity of finding this optimal 
MPS? It has been observed that DMRG converges expo- 
nentially fast to the ground state with a relaxation time 
proportional to the inverse of the gap A of the system 
|23j . For translational invariant critical systems, this gap 
seems to close only polynomially. As we have proven that 
D only have to scale polynomially too, the computational 
effort for finding ground states of 1-D quantum systems 
is polynomial (P) . Let us recapitulate for which systems 
and under which conditions this statement is true: 1) 
the a-entropy of blocks in the exact ground state grow 
at most logarithmically with the size of the block for 
some a < 1; 2) the gap of the system scales at most 
polynomially with the system size; 3) given a gap that 
obeys condition 2, there exists an efficient DMRG- like 
algorithm that converges to the global minimum. As the 
variational MPS approach [2£| is essentially an alternat- 
ing least squares method of solving a non-convex prob- 
lem, there is a priori no guarantee that it will converge 
to the global optimum |41|, although the occurrence of 
local minima seems to be very unlikely (2j|. Note that 
this e.g. implies that there would be no exponential gain 
in using a quantum computer in simulating 1-D spin sys- 
tems, as the steps needed in a quantum computer are 
also bounded by the inverse of the gap! 

Reduced density operators of translationally 
invariant states 

Now we turn to the study of how MPS approximate 
the local properties of translationally invariant states ip 
of N — > 00 systems. Let us denote by ip a translationally 
invariant state with N — 00 and by a — tr(p) the corre- 
sponding two-particle reduced density operator, where 
the trace is taken with respect to all particles except 
for two neighbors (say particle 1 and 2). Our goal is to 
show that there exists a Z3-dimcnsional MPS whose cor- 
responding reduced density operator, od, approximates 
well a. More specifically, we want to bound ||<r — 0\d||i 
from above. Note that any nearest neighbor correlation 
function, and in particular the energy density will be au- 
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tomatically bounded. 

One can represent the state ip |3^| as a MPS @ with 
all the A's equal and D,N = oo. Without loss of general- 
ity we can again choose the gauge condition (2) with A, A 
site independent. In fact, the elements An coincide with 
the eigenvalues of the reduced density operator corre- 
sponding to half a chain. As before, they will play a very 
important role in the bound that we derive; in particular, 
we will need e(D) = J2^Ld+i anc ^ we define 



We divide r and s in two parts, 



D 



D oo 



^0+Dn 
J 



7*= EVMEE 

\q=1 / 0=ln=l 



which roughly behaves like e(D), i.e. it becomes small 
when D is large. 

Lemma 3 Given ip and a as defined above, there exists 
a D -dimensional MPS such that 



^Tt\<j-<j d \ < 2(v/^D)+e(D))+ 7D 
Proof: One can easily show that 



(6) 



r i,i ,3,3 — r h* ,3,3 i _M ,3,3 
' — '0 ' ' 1 » 



The first parts (tq and sq) contain just the part of the 
sums with n = m = 0, whereas the second ones (fi and 
Si) contain the rest of the sums. 

The goal is to find upper bounds for 

Ar = J2 

As = Yl \^'' jd '\' 

At = j2 \rr' d ' j ' -sr' d ' j '\. 

We start with Ar. We can write 

d 

Ar= \Tr[AttAA i PA if A i 'i]+Ti[PAftAA i A if A3'i] 



r *,» ,3,3 = 



(i\<i'\p\j)\j') = Tr 



A^AA i A i ' A'"* 



(8) 



— ^ [A^ AA l ] a+Dn j j+Dm [A l A j ^](3+Dm,a+Dr, 
a,0,n,m 

Let us consider the Z?-dimensional matrix product den- 
sity operator (MPDO) 34] 

p D = Tr • • ■ B^°°) \h) . . . |U(ii| • ■ ■ (iool 

parameterized by 

p>i3 \" &i A*3 „ „ , 

c 'aa'\00 l ~ / , ^a+Dn.H+Dm^-g' +Dn.8' +Dm- ha , nd, a ,™ 



Here. 



P=X>)( 3 



x=D 



with \x) a unit vector in the computational basis. Using 
Cauchy-Schwartz inequality we have 



E \* 



A^ AA i PA % ' A^ 



< 



The following eigenvalue equations are easily verified: 

E B aa';0l3'^al3 = A a > pi 

where 



i,*',j;j'=l 



^ Tr^'U^A^'^"] 



1/2 



1/2 



d 2 [Tr[PAP]] 



1/2 



A a /3 — 5 a p 2^ 



The second term in Eq.© gives the same, so that we 
have 



Dk- 



k=0 



Ar < 2d 2 



The reduced density operator s" JJ = (i\(i'\pD\j)\j') is 
given by 



J2 K = 2d 2 ^/7(D). 

\ x=D 



a0nm 



+Dn,0+Dn 



A 1 ' 



(7) 



The term As is a bit more tricky. We write it as 
/ d 

A q — T, U ' T33* 

n,?n x,y i,i',j,j'= 1 



7 



where the prime in the sum means that n = m = is 
excluded, and 



T 11 

x,y,n,m 



A x ^ x f3,mA l x f3+Dn A l i 3A 



Dm.y ' 







Using Cauchy-Schwartz inequality we get 

As < EE E i<,.*.mi 2 = 

n,m x.y — 1 

d2 EE 

n,m /3 



Lfe A/3+Dfc 



— ^ 2 ^ ^ A^+Dn 
n=l 



< 2d 2 ^2\ x = 2d 2 e(D). 



=£> 



Now, it only remains to bound At. Again making use 
of the Cauchy-Schwartz inequality and the inequality 



1 



1 1 



1+xl+y 



< x + y 



we get 



D 



D oo 



a=l /3=ln=l y^P 

The bound states in the lemma is given by Ar + As + 
At. □ 

The usefulness of this lemma relies in the fact that the 
spectrum {A 1 } is typically decaying very fast for noncrit- 
ical systems, and hence the upper bound decays typically 
very fast with increasing D. Let us illustrate this on the 
hand of e.g. the XXZ Hamiltonian on a chain 

H = - Yl s ' S j + s ? S j + &S?S*. 

<i,j> 

For A < — 1, the spectrum of the halve-infinite chain can 
be calculated exactly [2^] . All eigenvalues are of the form 
c(z)z n where 

z = cxp (— 2arcosh(| A|)) 

00 

l/c(z) = J2a{n)z n 
n=0 

and where a(n) denotes the degeneracy of each eigenvalue 
and is given by the number of possible partitions of n into 
odd parts. The asymptotics of a(n) are known exactly 



course), we have solve the equation 

n(D) 

D = a(m) 



for n(D). Asymptotically, this becomes 

31/4 e ikV^ (\og(D)\ 2 

e(D), the sum of all eigenvalues smaller than the D 
largest ones, can now readily be calculated as 



e(D) = c(z) 



E 

m—n(D) 



a(m)z r ' 



exp 



3|logz| 



This last expression decays faster than any inverse power 
of D. A similar calculation holds for 7^, showing that 
the upper bound in lemma 3 decays faster than any poly- 
nomial in D. This is a very nice illustration for the su- 
perpolynomial accuracy of the infinite DMRG method 
for the XXZ-model. The other integrable models will 
have a very similar behavior, and it is expected that non- 
integrable models share the same features. Of course, the 
upper bounds derived here are not tight, and in practice 
the accuracy of a DMRG will be much better than de- 
rived here. 



CONCLUSION 

In conclusion, we highlighted the importance of the 
concept of local reduced density operators of transla- 
tional invariant systems, their approximation using ma- 
trix product states, and the connection with frustration 
and monogamy properties of entanglement. We quanti- 
fied how well matrix product states approximate exact 
ground states of 1-D quantum spin systems as a function 
of the entropy of blocks of spins, and showed that the 
complexity for representing ground states of 1-D systems 
(even critical ones) as MPS scales polynomially in the 
number of spins. 
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a(n) 



1 eTT 



4 3 1 /%3/4' 



To determine the power n(D) of the (D + l) th largest 
eigenvalue c(z)z n {D) (taking into account degeneracy off 
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